Additive Equivalence in Turbulent Drag Reduction by Flexible and Rodlike Polymers 
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We address the "Additive Equivalence" discovered by Virk and coworkers: drag reduction af- 
fected by flexible and rigid rodlike polymers added to turbulent wall-bounded flows is limited from 
above by a very similar Maximum Drag Reduction (MDR) asymptote. Considering the equations 
of motion of rodlike polymers in wall-bounded turbulent ensembles, we show that although the mi- 
croscopic mechanism of attaining the MDR is very different, the macroscopic theory is isomorphic, 
rationalizing the interesting experimental observations. 



I. INTRODUCTION 

Turbulent flows in a channel are conveniently discussed 
for fixed pressure gradients p' = —dp/dx where x, y and 
z are the lengthwise, wall-normal and spanwise directions 
respectively [1]. The length and width of the channel are 
usually taken much larger than the mid-channel height 
I/, making the latter a natural re-scaling length for the 
introduction of dimensionless (similarity) variables. Thus 
the Reynolds number Re., the normalized distance from 
the wall 2/+ and the normalized mean velocity V^{y^) 
(which is in the x direction with a dependence on y only) 
are defined by 



Re \ 



y 



yRe/L , 



(1) 



where is the kinematic viscosity. Drag reduction 
by polymers in a channel geometry is bounded by two 
asymptotes [2]. One is the von-Karman log-law of the 
wall for Newtonian fluids 



V+{y+) = n-^\ny+ ^ B , for j/+ > 30 



(2) 



While the log-law can be derived using several ap- 
proaches, the von-Karman constant Kk ~ 0.436 and in- 
tercept B « 6.13 are only known from experiments and 
simulations [1, 3]. The second asymptote is the Max- 
imum Drag Reduction (MDR) where the velocity field 
assumes another log-law of the form 

V+{y+) ^ —\n{eK^y+) for ?/+ > 10 . (3) 

This law, which had been discovered experimentally by 
Virk (and hence the notation k^) was derived theoret- 
ically for flexible polymers in [4-7]. The actual veloc- 
ity profile in the presence of polymers is bounded be- 
tween these two asymptotes: for sufficiently high values 
of Re and concentration of the polymer, the velocity 
profile in a channel is expected to follow the law (3). 
For finite Re , finite concentration and finite extension of 
the polymers one expects cross-overs back to a velocity 
profile parallel to the law (2), but with a larger mean 
velocity (i.e. with a large value of the intercept B). The 



position of the cross overs are not universal, and they are 
understood fairly well [7, 8]. 



In this paper we address the experimental finding that 
rigid rodlike polymers appear to exhibit a very similar 
MDR (3) as flexible polymers [9]. Since the bare equa- 
tions of motion of rodlike polymers differ quite signifl- 
cantly from those of flexible polymers, one needs to ex- 
amine the issue carefully to understand this similarity, 
which was termed by Virk "Additive Equivalence" . The 
aim of this paper is to understand this Additive Equiva- 
lence on the basis of the equations of motion. 



In Sec. II we address the theory of drag reduction by 
rigid polymers. In Sec. II A we consider the equations 
of motion of rigid polymers (or fibers) in the presence of 
strong wall-bounded turbulence. We explain the inter- 
esting differences between the interaction of flexible and 
rodlike polymers with turbulent fluctuations. In Sec. II B 
we discuss the difference in statistics between flexible and 
rodlike polymers near thermodynamic equilibrium. Sec- 
tion II C is devoted to the statistics of rodlike polymers 
in turbulent flows with a strong shear. In Sec. II D we 
address the important issue of how to evaluate the vari- 
ous cross-correlation functions between the polymer con- 
formation tensor and the turbulent fluctuations. These 
objects have a seminal role in the theory of drag reduc- 
tion by rodlike polymers. In the following Sees. HE and 
II F we demonstrate that in spite of the very signiflcant 
microscopic differences between flexible and rodlike poly- 
mers, the balance equations for momentum and energy 
have exactly the same form as the corresponding equa- 
tions for flexible polymers. The "additive equivalence" 
follows from this observation. In Sect. Ill we summa- 
rize the paper and discuss further the correspondence 
between drag reduction by flexible and rodlike polymers. 
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II. THEORY OF DRAG REDUCTION BY 
RODLIKE POLYMERS 

A. Basic equation of motion 

1. Hydrodynamic equations for the polymeric solutions 

The hydrodynamic equations for an incompressible 
fluid velocity U = U{t, r) in the presence of rodlike poly- 
mers have the form 



DU 
'Dt 

= V-U . 



vqAU - Vp + V • cr , 



Here the fluid density is set to unity {g 
the substantial derivative 



D 

Dt 



d_ 

dt 



(4a) 
(4b) 

1), D/Dt is 



(5) 



p = p{t, r) is the pressure field, vq is the kinematic vis- 
cosity of the carrier fluid and <r Uab is an extra stress 
tensor caused by the polymers. 

The calculation of the tensor cr for rigid rods is offered 
in the literature [10], subject to the assumptions that the 
rodlike polymers are mass-less and having no inertia. In 
other words, the rodlike polymers are assumed to be at 
all times in local rotational equilibrium with the velocity 
field. Thus the stress tensor does not have a contribu- 
tion from the rotational fluctuations against the fluid, but 
rather only from the velocity variations along the rodlike 
object. Such variations lead to "skin friction", and this 
is the only extra dissipative effect that is taken into ac- 
count. The result of these considerations is the following 
expression for the additional stress tensor: 

(Jab = Qi'pnanb (nirijSij) , rodlike polymers , (6) 

where Vp is the polymeric contribution to the viscosity 
at vanishingly small and time-independent shear; Up in- 
creases linearly with the polymer concentration, making 
it an appropriate measure for the polymer's concentra- 
tion. The other quantities in Eq. (6) are the velocity 
gradient tensor 



S = {VUy => Sab = dUa/dXb , 



(7) 



and n = n{t,r) is a unit {n ■ n = 1) director field that 

describe the polymer's orientation. Notice, that for flex- 
ible polymers the equation for aab is completely different 
from Eq. (6): 



<Jab = i'pTp naUb , flexible polymers . 



(8) 



Here 7p is the polymeric relaxation frequency. The differ- 
ence between Eqs. (6) and (8) for aab for the rodlikc and 
flexible polymers reflects their very different microscopic 
dynamics. For the flexible polymers the main source of 
interaction with the turbulent fluctuations is the stretch- 
ing of the polymers by the shear. This is how energy 



is taken from the turbulent field, introducing an addi- 
tional channel of dissipation without necessarily increas- 
ing the local gradient. In the rigid case the dissipation 
is only taken as the skin friction along the rodlike poly- 
mers. Having in mind all these differences it becomes 
even more astonishing that the macroscopic equations for 
the mechanical momentum and kinetic energy balances 
are isomorphic for the rodlike and flexible polymers, as 
is demonstrated below. 



2. Equation for the orientation field of dilute rods solution 



The equation for n{r,t) has the form 

— = (5 - nn) • 5 • n => = {6ai - naUi) S^rij . 

(9) 

This equation conserves the unit length of the director 

n. The theory is written more compactly in terms of the 
(normalized) conformation tensor 



72. = nn => TZab = riaUb 



(10) 



The equation of motion of this object follows from 
Eq. (9): 



DUab 

Dt 



= SaiR-ib + SbiR-ia - 2 TZabi^ijlZij) . (11) 



In terms of TZ wo can rewrite cr, Eq. (6), as 
a = 6upnTr{n- S} ^ Uab = &yp1lab{1lijSij) . (12) 



B. The statistics of rodlike polymers near 
thermodynamic equilibrium 

In order to understand the interaction of rodlike poly- 
mers with wall bounded turbulent fluctuations, we need 
to start along the well charted road of dynamics near 
equilibrium, and examine the solution of these equation 
for strong shears. 



1. Equations of motion near thermodynamic equilibrium 

Consider the shear S{r,t) in Eq. (9) as independent 
of space and time, and replace S{r,t) =J> S. Near equi- 
librium rodlike polymer (or fibers) experience rotational 
disorder due to local thermal velocity fiuctuations that 
can be c;onsidorod as a Brownian motion in the space of 
angles. As usual, the effect of thermal fiuctuations can 
be mimicked by adding to the RHS of Eq. (9) a Langcvin 
random force f{t). When the rigid polymer is symmet- 
ric the hydrodynamics can depend only on the dyadic 
product nn. By pre- and post-multiplying the equation 
above by 715, and taking average over the Brownian fiuc- 
tuations, we have the transport equation for the second 
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moment Rab = ^^a^^6, where and ovcrbar indicates an 
average over the polymer configurations; 

— = Sai'R-ib+Sbi'R-ia-'iSijTZijTZab-QjB (^ab ^ 

(13) 

where 7b is the Brownian rotational frequency, propor- 
tional to the temperature. The derivation of this equa- 
tion in the literature employs the closure 



T^ijl^ab ^ab — '^ij'^ab ^ab ; 

which is rationalized in [10] 



(14) 



where the angular brackets denote an average over the 
turbulent fluctuations. For well developed turbulence it 
is expedient to use the Reynolds decomposition in which 
the velocity field U{r) is written as a sum of its average 
(over time) and a fluctuating part: 

U{r, t) = V{y) + u{r, t) , V{y) = {U{r, t)) . (22) 

For a channel of large length and width all the aver- 
ages, and in particular V{y) ^(2/)j are functions of 
y only. Correspondingly the shear is written as a sum 
of the mean shear Sab and fluctuating shear s(r, t) with 
zero mean: 



2. Solution for a simple shear flow 
For a simple shear the velocity gradient satisfies 

Sab = SSaxSby , (15) 

Eq. (13) in the stationary, space-homogeneous case turns 
into: 

-DCr {SaxT^yb + hxT^ya — '^T^ab'^xy) = 2(37?.a6 — Sab) ■ 

(16) 

Here we introduces the Deborah number for the rod-like 
polymers: 



(17) 



Equation (16) was solved in the limit Dbj. 1 [10] with 
the final results 



Tlx 



27? ^ 

^ '^xy 



1 > Tlxy 



2Der 



1/3 



21/3 



De 



2/3 ' 



(18) 



It is interesting to compare the statistics of rodlike and 
flexible polymers in strong laminar shears. For example 
(cf. [6] and references therein). 



Tlx 



xy 



Del'^Tlyy , rodhke , 



Dei Tly 



flexible 



(19) 
(20) 



Here the Deborah number for flexible polymers Dei is de- 
flned with the flexible polymer relaxation frequency 7p, 
Dei = S/'jp. The different dependence on Deborah num- 
ber stems from the very different microscopic dynamics 
that leads to different expressions for the polymeric stress 
tensor aab for the rodlike and flexible polymers. 



C. Statistics of rigid rods in turbulent flow with 
strong shear 

In the presence of turbulence the fluctuations arc no 
longer thermal and the statistical description of the poly- 
mer orientation is accomplished with the mean values 



Rab = {Tlab) 



(21) 



S = {S{r,t)), S{r,t) = S + s{r,t), {s{r,t)) = 0, 

(23) 

With these notations and after averaging over the statis- 
tics of turbulence, Eq. (11) takes the form 



DTZab 

Dt 



Aab - B. 



ab ) 



(24) 



where Aab contains all the terms in which the mean shear 
appears explicitly, and Bab contains only the fluctuating 
part of the shear: 

-^ab = SaiRib + ^biRia ~ 2 Sij {TZabTlij) , 

Bab = {Tlab{SijTZij)) - {SaiTZib) - {SbiTlia) -(25) 

In a steady state in channel geometry the LHS of Eq. (24) 
contains only one non- vanishing terms, which is (it- V7?.). 
This term is responsible for the turbulent part of the 
physical flux of TZ. In all our derivation below we will 
assume that such terms are small compared to all 'local' 
contributions to the balance equations for the momentum 
and energy. This assumption will have to be tested a- 
posteriori. 

With this in mind Eq. (24) reads simply Aab = Bab- 
Thin will allow us to estimate the crucial correlations 
functions that appear in the theory of drag reduction 
below. The LHS of this relation can be made explicit in 
channel geometry; using Eq. (15) we flnd 



Aab = S {6axRyb + ^bxRya - 2 {TZabT^xy)) ■ (26) 

Using definition (10) and constraint \n\ = 1 this Eq. can 
be rewritten in components as: 



A-XX 2S (J^xyiT^yy -\- 7^22 )) 

A-yy — 2 (5' (T^xyT^yy) , 

Azz = —2 5' (TZxyTi-zz) , 

Axy S (7?-yy(l 2 7^3,3,)) . 



(27a) 
(27b) 
(27c) 
(27d) 



In writing down expressions for Bab we will make ex- 
plicit use of the expected solution for the conformation 
tensor in the case of large mean shear, S'^ 3> (s^). In 
such flows we expect a strong alignment of the rod-like 
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polymers along the streamwise direction x. The direc- 
tor components Uy and are then much smaller than 
Tlx ~ 1. For large shear we can expand according to 



1 



1 



(28) 



Wc note that for nix = 1 (when the shear is actually infi- 
nite) the object Bab vanishes since (s) = 0. We therefore 
represent Bab in order on Uy Uz, keeping up to second 
order. We will show below that it is important to keep 
terms of second order since some of them have the same 
magnitude as terms which are formally of first order in 
the smallness. All terms of third order are smaller in 
magnitude than the terms that we keep. The first two 
orders read: 



yy 

3oy 



2 {^xy'^xy H" ^xz^xz) •> (29a) 

-2 {s^yTl^y) , (29b) 

-2 (s,,7^,,) , (29c) 

{'^xyi.^xx ^yy) "^xz^yz) •> (29d.) 



^xx — ^ i'^yyi'^'^yy '^zz) '^'T^zzi^^zz ^xa;) 



i.yyy^.Dyy 



)> 



(30a) 



B(2) 

yy 



(J^yy i^Sxx ^S^^y ) 2 T^xZ-^XzY) ' (30b) 

(j^yyi^Sxy + ^Syx ) + T^zzSyz (30c) 



Equations (27), (29) and (30) serve as basis for further 
analysis. 



D. Closures and orders of magnitude 

1. Statistical objects of interest 

A theory of turbulent channel flows of Newtonian fluids 
can be constructed in terms of the mean shear S{y), the 
Reynolds stress W{y) and the kinetic energy K{y); these 
are defined respectively as 

S{y) = dV{y)/dy , W{y) = -{u^Uy) , K{y) = {\u\^)/2. 

(31) 

In the rodlikc polymer case the additional stress tensor 
aij and its various correlation functions needs to be con- 
sidered as well. For that purpose we turn now to the 
analysis of the necessary statistical objects. 

First note that the expansion (28) allows us to express 
all products TZabT^cd = naUbncTid in terms that are linear 
in 72., up to third order terms in Uy Hz- For example, 



7? 7? 

'^xy '^xz 



1 - 2{ny 



7^: 



xy 



■yy 5 (32) 
etc. 



Actually we have used these estimates in the derivation 
of Eqs. '(29) and (30). 



As a further preparation for the theory below we ana- 
lyze various statistical objects in the turbulent environ- 
ment and estimate their magnitudes. Based on experi- 
mental observations and DNS data we assume that statis- 
tics of turbulent fluctuations does not deviate too much 
from isotropy. Explicitly, 



i^xx) 



yy/ 



{Slz) 



i^xy) 



(33) 



Here and below the notation means "the same order of 
magnitude (i.e. correct to leading order up to coefficients 
of the order of unity)" . 

Second, consider correlation functions of turbulent 
fluctuation of the shear, s^, i.e. {sijSi~e)- At distance 
y from the wall the correlation functions is dominated by 
eddies of size y. Thus 



{sijSke) ^ K{y)/y^ , for aU ij, kl 



(34) 



Third, consider a cross-correlation functions of the ten- 
sors TZ with the turbulent shear s. For this goal take the 
leading terms in the LHS and the RHS of xy Eq. (27d) 
and (29d): 



\SR 



K 



■yy\ 



r(i)| < B — — 

^xy I ~ -"-xj/ ) 



y 



(35) 



where on the RHS we have used the Cauchy-Schwartz in- 
equality and Eq. (34). The notation < means "~ in the 
sense defined above, or having smaller order of magni- 
tude". Below, (cf. Sec. HF) we show that this estimate 
is saturated. Therefore expecting that Ryy ~ R^y we 
have 



SR 



xy 



K 

y 



(36) 



This estimate can be recast into an intuitive form which 
is 



Rxy — (^x^y) 



(37) 



This is in direct accord with the understanding that the 
degree of deviation from perfect alignment of the rod- 
like polymers {Rxx = 1) is proportional to the turbulent 
fluctuations relative to the mean shear. 

Taking now the leading terms in yy Eq. (27b) and (29b) 
we have 



S IjZyyTZxy) — 2^^^^ ^ {SxyTZxy) 



(38) 



Estimating {TZyyTZxy) as ~ RyyRxy, and using Eq. (36), 
we have 



(39) 



{.^xy'^xy^ ^ -^yy * 



Notice then the Cauchy-Schwartz inequality for the same 
correlation, 



{Sxy'^xy) ^ ^x 



(40) 
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gives much higher upper bound, than the real estimate 
Eq. (39). This shows how important is the use of the 
equations of motion in estimating various correlation 
functions; sometime the direct Cauchy-Schwartz estimate 
saturates, and sometime it is a gross overestimate. The 
reason for why different correlations have different order 
of magnitude can be traced back to Eqs. (27) which in- 
dicate that the diagonal components of A are cubic in 
the small parameter Uy ^ riz while A^y is quadratic. 

Notice that the correlator {TZijSij) has the contribu- 
tions of the type of {sxy'R'xy) and {syyTZyy). Both of 
them have the same estimate RyyVK/y. Therefore we 
can write 

{TZijSij) Ryy^ . (41) 



E. The momentum-balance equation 

At this point we can apply our estimates in the context 
of the balance equations for the mechanical momentum 
and the energy. We begin with the former, which is exact. 
It reads 



uoS + W + {a^y) = p'{L - y) 



(42) 



where p' = —dp/dx. Near the wall y L and the RHS 
of this equation is approximated as p' L, a constant pro- 
duction of momentum due to the pressure gradient. On 
the LHS we have the Reynolds stress which is the "tur- 
bulent" momentum flux to the wall, in addition to the 
viscous and the polymeric contributions to the momen- 
tum flux. 

Using Eq. (6) and Reynolds decomposition (23) we 
compute 

(o-xy) = GiypiUxyTZijSij) = 6 Up[S (Uly) + {TlxyTZijSij}] . 

(43) 

With Eq. (32) the first term in the RHS of this equation 
can be estimated as follows: 

6 UpS (nly) = ClUpRyyS , ci ~ 6 . (44) 

On the other hand, using the estimate (41) one sees that 
the second term in the RHS of Eq. (43) is negligible. 

Finally we can present the momentum balance equa- 
tion in the form 



voS + ciiypRyyS + W= p'L 



(45) 



Another way of writing this result is in the form of an 
effective viscosity. 



u{y)S + W=p'L 



(46) 



where the effect of the rodlike polymers is included by 
the effective viscosity z^(y): 



Hy) = 1^0 + ciUpRy 



(47) 



F. Turbulent Energy Balance Equation 

In considering the balance of energy in a channel flow 

it pays to separate the spatial directions, since wo can 
learn separate bits of information from each such equa- 
tion. Introduce the partial kinetic energy density 

K,{y) = ^{ul), K{y)=Kx + Ky + K, , (48) 

and consider the partial energy balance of Ka{y) 
dKa{y) 



dt 



+ Ra+ sT + £P = Wiy)Siy)6ax • (49) 



The total density of the kinetic energy at given distance 
y from the wall is: 



9K{y) , ,dis , n_ 



dt 



+ s'''- + sP = W{y)S{y), 



(50) 



The various symbols in the last two equations are ex- 
plained as follows: The RHS of these equations de- 
scribes the energy flux from the mean flow to turbulent 
fluctuations due to the correlation between stream-wise 
and cross-stream components of the turbulent velocity, 
known as the Reynolds stress W, see Eq. (31). Remark- 
ably, in channel geometry this flux exists only in the equa- 
tion for the streamwise velocity fluctuations, K^- 

The term Ra (y) is known as the "Return to isotropy" 
[11], and it vanishes for isotropic turbulence in which 
Kx = Ky = Kz- Otherwise it redistributes partial kinetic 
energy between different vectorial components and does 
not contribute to the total balance (50). A simple model 
for this term [11] is 



Ra 



K 

y 



3 J 



(51) 



As usual the local "outer scale of turbulence" was esti- 
mated as the distance to the wall y. The order of magni- 
tude estimate (33) is in accord with the role of this term 
in returning to local isotropy. 
The term 



^a" - ^0 X] (^?a) ' no sum over a , 



(52) 



on the LHS of Eq. (49) is the rate of the viscous dis- 
sipation, proportional to the kinematic viscosity of the 
carrier fluid uq. Lastly, the polymer contribution to the 
energy balance, denoted as s^, can be exactly computed 

as 

£a = {<^ajSaj) = 6 f p {SaiHai {SUxy + SjkTljk)) ■ (53) 

Notice that Eqs. (49) and (50) are written "in the local 
approximation" , in which the energy flux in the physical 
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space is neglected. This is consistent with neglecting the 
term {u ■ V7^) in our discussion after Eq. (25). A jus- 
tification of this approximation in the problem of drag 
reduction by polymers is found in [4, 5]. 

Using the expansion (28) we can rewrite the equation 
for the dissipation rate eP as a series in the small param- 
eter Uy ^ n^: 

~ 6fp ^{Sa;3; + \T^xy{^xy "I" ^yx) "I" ^xzi^xy ~^ ^j/x)] i 

+ [syyilyy + SyzHyz + more] 3 + ' ' ' } (54) 

><{sn^y + s^^ + [---]^ + [---]^ + ...]) . 

In the square brackets [•••]! we displayed all the terms 
which are linear in the small parameter. In the square 
brackets [• • • ]2 we show two quadratic terms, and there 
are more of them as indicated. The symbol + . . . stands 
for "higher order terms". In the second curly brackets 
the square brackets [■ and [■ • • ]2 are identical to the 
corresponding terms in the first curly brackets. The two 
leading terms in Eq. (54) are proportional (s^^) and 
S (TZxySxx)- Using Eqs. (34), (29d) and (35) one sees 
that both leading terms have the same order of magni- 
tude, K{y)/y'^. In fact we will argue that these two 
terms must cancel each other, up to terms of the a higher 
order contributions --^ RyyK{y)/y'^. 

To see that such a cancellation must exist, we note 
that near the MDR we expect the polymer contribution 
to the dissipation to balance the production term WS (cf. 
Eq. (58) below). We will show later that this production 
term is y independent (where y is the distance from the 
wall). On the other hand we will show that K(y) is linear 
in y. making K/y^ very large near the wall. Therefore 
WS cannot be balanced by K/y"^. To avoid using the 
final results at this stage, and nevertheless to see that a 
cancellation must exist, we can focus just on the station- 
ary balance Eq. (49) for the y component, in which the 
RHS is zero. Notice also that near the MDR the poly- 
mer contribution to the energy balance dominates over 
the viscous dissipation and the nonlinear energy fiux from 
large to small scales [4-6]. The latter was evaluated in [4] 
as K^l"^ jy which is exactly the evaluation of the return 
to isotropy term. The conclusion is that near the MDR 

= {(^ajSaj) =&i^p {syiTZyi {STZxy + SjkR-jk)) « . 

(55) 

This expression can be again arranged in orders of mag- 
nitude, similarly to Eq. (54): 

£^ = ^{ [T^-xySj/x] yi + \^yyT^yy + SyzTi-yz] y2 } 

x{STl^y + s^^+[---]^ + [---]^ + ...}) . 

Here the second curly brackets is the same as the second 
curly brackets in Eq. (54). In Eq. (56) we find again the 
same two large contributions, i.e. TZxySyxiSTZxy + Sxx}+ 
1.0. t, where the lower order terms have at least one small 
factor. Thus these terms must cancel each other, which 
means that STZxv cancels s^x inside correlation functions. 



But these are exactly the terms that appear in the sum 
of the two large terms in Eq. (54), and we are therefore 
justified in neglecting them, allowing the other terms to 
share the burden of balancing WS. We note that this 
conclusion is justifying a-posteriori the statement after 
Eq. (35) that the Cauchy-Schwartz inequality is satu- 
rated for the case considered. 

An additional way to see that the cancellation must 
take place is to examine again the momentum balance 
equation (45). As discussed below, the MDR is obtained 
when the Reynolds stress term W is negligible compared 
to the polymer contribution ciVpRyyS. But when this 
happens it means that 

WS < ciUpRyyS^ « ciVpK/y'^ . (56) 

Evidently this means that also in the energy balance 
equations WS would be overwhelmed by terms of the 
order of K/y'^ which therefore must cancel against each 
other. 

Using our order of magnitude estimates for the remain- 
ing terms we can therefore conclude that 

s^{y)^C2ypRyy{y)K(y)/y\ (57) 

where C2 is another parameter of the order of unity. Re- 
turning to the balance equation for the energy, we recall 
that we cannot calculate e'*"'(y) exactly, but we can esti- 
mate it rather well at a point y away from the wall. When 
viscous effects are dominant, this term is estimated as 
i^{a/y)'^K{y) (the velocity is then rather smooth, the gra- 
dient exists and can be estimated by the typical velocity 
at y over the distance from the wall). Here a is a con- 
stant of the order of unity. When the Reynolds number 
is large, the viscous dissipation is the same as the turbu- 
lent energy fiux down the scales, which can be estimated 
as K{y)/T{y) where T{y) is the typical eddy turn over 
time at y. The latter is estimated as yjh^J K{y) where b 
is another constant of the order of unity. Together with 
Eq. (57) we can thus write the energy balance equation 
at point y as 

K(v) K^/'^iv) Kiv) 
av^^ + b ^ + c2iypRyy{y)^ = W{y)S{y) . 

We recognize the important result that the effective vis- 
cosity induced by the rodlike polymers in both the mo- 
mentum and the energy balance equation is proportional 
to Ryy. These balance equations are identical in form to 
those found for flexible polymers [4]; this is an impor- 
tant step in understanding the "Additive Equivalence" 
discovered by Virk. 

To complete the derivation one adds to the balance 
equation the relation between K(jj) and W{y) which in 
the elastic layer are expected to be proportional to each 
other, 

K = clW . (59) 

It should be stressed that rigorously one can establish 
this relation only as an inequality with Cv < 1, and its 
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use as an cqiiality (which is common to the derivation 
of the Newtonian log-law as well as to the derivation of 
the MDR in the flexible polymer case) rests on exper- 
imental and simulational confirmation. Near the MDR 
the terms representing the effect of the polymers in Eqs. 
(45) and (58) are dominant, and one estimates from the 
momentum equation Ryy{y) oc 1/S{y). Using this in Eq. 
(58) together with Eq. (59) one ends up with the pre- 
diction that S{y) oc leading to a logarithmic law 
for the mean velocity. Repeating the derivation of [4] in 
wall units one ends up with the MDR Eq. (3), with the 
identification 

= Cv/cn2/+ . (60) 

In this equation Cn and are constants that appear in 
the Newtonian theory, and cannot change from flexible 
to rodlike polymers. The existence of drag reduction is 
guaranteed, since Kv was shown to be larger than its 
Newtonian counterpart Kk [4]. The actual value of the 
slope at the MDR logarithmic law depends nonetheless 
on the numerical value of c^,. Thus the prediction of 
the theory is that if Cy is about the same in rodlike and 
flexible polymers, than the slope of the MDR should be 
about the same. 



III. CONCLUSIONS 

We have presented a scenario to rationalize the "ad- 
ditive equivalence" discovered by Virk and coworkers. 
The main conclusion of this paper is that although, on 
the face of it, the dynamics of flexible and rodlike poly- 
mers appear difl'erent, with flexible polymers being able 
to "stretch" and "store" energy (something that many re- 
searchers thought is central to drag reduction) , the effec- 
tive Reynolds balance equations for momentum and en- 
ergy are isomorphic. Accordingly, the MDR is expected 
to be the same as long as Cy of Eq. (59) is the same. 

One should stress however that Cv is not expected to 
be a universal number. It appears to us unlikely that Cy 
remains exactly the same in rodlike and flexible polymers. 
It appears more likely that these respective numbers are 
of the same order, giving the impression that the MDR 
asymptotes are the same. In addition, one expects that 
the cross-over from the MDR to the Newtonian plug, 
which is non- universal even in flexible polymers [7, 8], 
may show significant differences between rodlike and flex- 
ible polymers. Indeed, in friction coordinates drag reduc- 
tion by rodlike polymers appears as an upper bound on 
the drag reduction by flexible polymers [9] . According to 
the theory of [8] flexible polymers reach their maximal 
drag reduction when fully stretched, being then as effec- 
tive as rodlike polymers. This is one way of rationalizing 
the findings of [9]. 

To make the difference between the flexible and rodlike 
polymers sharper, we note the different y dependence of 
K{y) and W{y) in the two cases. In the flexible case one 
had a threshold condition for the onset of drag reduction 



in terms of the Deborah number, stating that the typical 
time scale for turbulent fluctuations, y / y^K (y) is of the 
order of the polymer relaxation time Tp. This immedi- 
ately leads to the estimate K{y) ~ y"^, and the same for 
the Reynolds stress. In the present case we have esti- 
mated SRxy ~ ^/K/y, and with R'^y Ryy ?/ we get 
K{y) ~ y, and due to Eq. (59) we can write 

K{y) ^ W{y) ^ y , for rigid rodlike polymers; 
K{y) W{y) y^ , for flexible polymers. (61) 

We note that this last statement is a sharp prediction for 
an important difference between the two drag reducing 
universality classes, a difference that is not at all in con- 
tradiction with the "additive equivalence" . A-posteriori 
we can also see why the terms of the order of K{y)/y^ 
in Eq. (54) must have cancelled, being divergent as 1/y 
against a y independent energy input W{y)S{y). We 
hope that this prediction would be put to experimental 
or simulational test. 

An additional important difference between rigid and 
flexible polymers is that in the latter case an important 
condition for attaining the MDR was Rxx ^ Ryy ■ In the 
present case we need to guarantee that ciVpRyy ^ z/q 
in order to enable the polymer terms to overwhelm the 
Newtonian terms in the balance equations. This condi- 
tion means however that the concentration of the rigid 
polymer should be large enough before the MDR is ob- 
tained. In the flexible polymer case one could reach the 
MDR conditions even for small concentrations as long 
as the Re is large enough and the Deborah number is 
large, leading to Rxx ^ Ryy [7]- This difference leads 
to the observed experimental behavior, where for flexible 
polymers the MDR is reached even for small concentra- 
tions and then a cross over back to the Newtonian plug is 
found, whereas in rigid polymers that MDR is obtained 
gradually as the concentration increases, and see the fig- 
ures in [9] for comparison. 

Finally, it is interesting to note that our order of mag- 
nitude estimates of Rij could be read directly from Eqs. 
(18) by replacing the thermal mean values TZij with tur- 
bulent means Rij = (R-ij) and simply identifying the 
Brownian frequency 7b in the definition (17) of the Brow- 
nian Deborah number De^ with the characteristic turbu- 
lent frequency 7turb = ^KRyy/y. Put for example in the 
equation for Ryy we get 

Ryy « 2''Ky'R%y[S{y)y]y' . (62) 

Simplifying, this equation reads S^Ryy ^ K{y)/y'^ which 
is nothing but the square of Eq.(36). All the other orders 
of magnitude derived in Sect. II D follow as easily with 
this identification. We believe that this is another way 
to argue that our estimates are physically sensible and 
that we have captured the essence of drag reduction by 
rodlike polymers and the nature of the observation of the 
"additive equivalence" . 
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